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Abstract We present a comprehensive review of MHD wave behaviour in the neigh- 
bourhood of coronal null points: locations where the magnetic field, and hence the local 
Alfven speed, is zero. The behaviour of all three MHD wave modes, i.e. the Alfven wave 
, and the fast and slow magnetoacoustic waves, has been investigated in the neighbour- 

ed . hood of 2D, 2.5D and (to a certain extent) 3D magnetic null points, for a variety of 

assumptions, configurations and geometries. In general, it is found that the fast mag- 
. netoacoustic wave behaviour is dictated by the Alfven-speed profile. In a /3 = plasma, 

' the fast wave is focused towards the null point by a refraction effect and all the wave 

energy, and thus current density, accumulates close to the null point. Thus, null points 
will be locations for preferential heating by fast waves. Independently, the Alfven wave 
is found to propagate along magnetic fieldlines and is confined to the fieldlines it is 
O ' generated on. As the wave approaches the null point, it spreads out due to the diverging 

fieldlines. Eventually, the Alfven wave accumulates along the separatrices (in 2D) or 
^ . along the spine or fan-plane (in 3D). Hence, Alfven wave energy will be preferentially 

' dissipated at these locations. It is clear that the magnetic field plays a fundamental 

role in the propagation and properties of MHD waves in the neighbourhood of coronal 
null points. This topic is a fundamental plasma process and results so far have also 
^ , lead to critical insights into reconnection, mode-coupling, quasi-periodic pulsations and 

' phase-mixing. 

Keywords Magnetic fields. Coronal • Magnetic fields, models • Waves, Magnetohy- 
. drodynamic • Waves, Propagation • Magnetohydrodynamics 
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' 1 Introduction 



X 



Magne tohyd rodynamic (MHD) wave motions (e.g. Roberts l2004l : Nakariakov & Ver- 
wicht e I2OO5I : De Moortel boosl ) are ubiquitous throughout the solar corona (Tomczyk 
et al. [20071). Several different types of MHD wave motions have now been observed by 
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various solar inst rumen ts: slow magnetoacousticwaves have been seen in SO HO data 
(e. g. Ofm an et aZ. 1 19971: D eForest & Gurman [l99i : Berglimans fc Clet te [l999l : Kliem et 
a/. l2002l : Wang et aZ. 120021 ') and TRACE data (De Moortel et a/ I2OO0I). Fas t magnetoa- 
coustic waves have been seen with TRACE (Aschwanden et a;. ll999l.l2o63 : Nakariakov 



20081'). Non-thermal 



et al. I1999I : 'Wang & Solanki liooi ') and ffinorfe (Ofman & 'Wang 
line n arrowing / broad ening due to Alfven wa ves ha s been reported by Harris on et al. 
(2003) / Erdelyi et al. (|l998l ). Banerjee et al. (|l998l 'l and O'Shea et a/. (20031 '). Alfven 



wa ves h ave possibly been observed in the cor ona (O kamoto et al. 200^ . Tomczvk et 
al. l2007l ) and chromosphere (De Pontieu et al. l2007l : Jess et al. 2OO9I ). although these 



cla ims a re subject to intense discussion (Erdelyi & Fedun 120071 : Van Doorsselaere et 
al. I2OO8I V 

It is clear that the coronal magnetic field plays a fundamental role in the propaga- 
tion and properties of MHD waves, and to begin to understand this inhomogeneous, 
magnetised environment, it is useful to look at the topology (structure) of the magnetic 
field itself. Potential-field extrapolations of the coronal magnetic field can be made from 
photospheric magnetograms, and such extrapolations show the existence of important 
features of the topology: null points - locations where the magnetic field, and hence the 
Alfven speed, is zero, and separatrices - topological features that separate regions of 
differe nt magnetic flux connectivity. A comprehensive review can be found in Longcope 
( 20051 ). 

This paper will provide a comprehensive literature review of the nature of MHD 
wave propagation in the neighbourhood of coronal null points. This topic exists at 
the overlap of two important areas of solar physics: MHD wave and magnetic null- 
point theories. A brief introduction to both of these areas is provided in i)1.2l and 
i)1.3l including a mathematical description of magnetic null points. ij2] reviews the 
early work that considers a 2D null point in a cylindrically symmetric geometry and 
describes the system in terms of normal modes. ^ reviews work performed in a 2D 
cartesian geometry that focuses on externally driven perturbations, and iJS. 61 describes 
the extension of these investigations into the nonlinear regime. SJ?] details the effects 
of threading the 2D X-point with an orthogonal weak-guiding field. iJ5] details the 
behaviour of MHD wave propagation in the neighbourhood of 3D null points, and the 
conclusions and summary are given in ij6l 



1.1 MHD Equations 

The viscous, resistive, compressible MHD equations utilised in this paper are: 



dt 



+ (v- V)v 



= -Vp + - (V X B) X B + !/V ■ TT 



^ = V X (v X B) + ryV^B 
at 



dp 
dt 
de 



+ V • (/9V) = 



dt 



+ (v- V)e 



(1) 



where p is the mass density, v is the plasma velocity, B the magnetic induction (usually 
called the magnetic field), p is the plasma pressure, /i = 47r x 10~^Hm~^ is the magnetic 
permeability, v is the coefficient of classical viscosity, tt^ 



ij — SijV ■ V is the stress 
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tensor, Sij — [dvi/dxj + dvj /dxij /2 is the rate-of-strain tensor, a is the electrical 
conductivity, 77 — 1/i-iu is the magnetic diffusivity, e — p/p{'y—l) is the specific 
internal energy density, where 7 = 5/3 is the ratio of specific heats and j = V x B//i 
is the electric current density, v and 77 are assumed to be constants. 

Note that the classical viscous term used in equations ([TJ is in fact not the most 
appropriate for the solar corona since, in the presence of strong magnetic fields, the 
viscosity takes the f orm of a no n-isotr opic tensor. However, only the papers of Craig 
& Litvinenko (|2007l l and Craig (|2008l ) mentioned in this review will invoke the non- 
isotropic viscous tensor and so, for brevity, we do not provide a full description here. 
The m athematical details of the non-isotropic v iscous tensor can be f ound in Braginskii 
(|l965l ) and, f or example. Van der Linden et al. l|l998l ). Ofman et al. (|l994l ) and Erdelyi 
& Goossens (|l994l : Il995i '). 



1.2 MHD waves 

A wave is a disturbance that propagates through space and time, usually with the 
transference of energy. Such a disturbance, either continuous or transient, propagates 
by virtue of the elastic nature of the medium. In MHD, the magnetic tension provides 
an elastic restoring force, such that we would expect waves to propagate along uniform 
magnetic fieldlines with a characteristic speed: 

|B| 

VA 



where is called the Alfven speed. Transverse waves travelling at this speed along 
magnetic fieldlines are called Alfven waves. 

If we consider a compressible medium, then we can define the sound speed as: 



Cs = 

When assuming a compressible medium, the Alfven wave still remains, but the sound 
and Alfven speed can now couple together to give magnetoacoustic waves. Two com- 
binations arise: the higher frequency mode is known as the fast magnetoacoustic wave 
and the lower frequency wave is known as the slow magnetoacoustic wave. These three 
wave types, the Alfven wave and the fast and slow magnetoacoustic waves, make up 
the three MHD waves considered in this review paper. 

The fundamental properties and nature of linear MHD waves in uniform magnetic 
fields have been reported in de tail by several authors, for example in an unbounded 
homogeneous medium (Cowling [l97y)i a-nd in a bounded inhomogene ous sl ab / cylin- 
drical den sity pr ofile embedded in a uniform magn etic fi eld (Roberts [l98ll ') / (Ed win 
& Roberts ITgia : Callv ll986l : Roberts & Nakariakov boosh . 

Finally, in MHD it is important to consider the ratio of magnetic pressure to thermal 
pressure. This ratio is called the plasma—/? and is given by: 

|B|^ 

The properties of the fast and slow magnetoacoustic waves have a strong dependence 
on the magnitude of the plasma—/?, namely because it is directly proportional to the 
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square of the ratio of the sound speed to the Alfven speed. Thus, in a regime where 
/? ^ 1, magnetic pressure and magnetic tension dominate tlie propagation and vice 
versa. Table 1 lists the main properties of the three wave types depending upon their 
environment. Note that the Alfven wave behaviour is independent of the plasma-/?, as 
it is a purely magnetic wave (in the linear regime). 

The plasma—/? parameter varies greatly with height across the layers of the solar 
atmosphere (see Garv [ioOll for well-constrained values). However, magnetic pressure 
generally dominates thermal pressure in the solar corona, and thus it is usual to assume 
plasma—/? <^ 1 when modelling a coronal environment. Hence, when we talk about 
coronal null points, we are talking about null points in a low— or zero—/? environment, 
although there are some caveats to this ( i]3.5p . 





plasma-/3 > 1 (high-/3) 


plasma—/? <C 1 (low-/?) 


Alfven wave 


Transverse wave propagating at speed va 


Fast MA wave 


Behaves like sound wave 
(speed Cs) 


Propagates roughly isotropically 
Propagates across magnetic fieldlines 
(speed va) 


Slow MA wave 


Guided along B 
(speed Va) 


Guided along B 
Longitudinal wave propagating 
at speed Cs 



Table 1 Properties of MHD waves depending on the plasma—/?. 



1.3 Magnetic Topology 

The magnetic field plays an essential role in understanding the myriad of phenomena 
in the solar corona. A realistic magnetic field can have many different components, 
and we can use topology nomenclature to reduce a complicated set of fieldlines to 
something more understandable. In 2D, a general magnetic configuration contains sep- 
aratrix curves (separatrices) which split the magnetic plane into topologically distinct 
regions, in the sense that within a specific region all the fieldlines start at a particular 
source and end at a particular sink. There is a second important topological aspect: 
magnetic null points (or neutral points) are single-point locations where the magnetic 
field vanishes (B = 0). There are two types of magnetic null point: X-type null points, 
commonly called X-points, which occur at the intersection of separatrix curves, and 
0-type null points, or 0-points, located at the center of magnetic islands. Magnetic 
topologies that contain null points are common in the presence of multiple magnetic 
sources. 

A magnetic fieldline that joins two null points (itself a special type of separatrix) 
is called a separator. Thus, instead of showing all the magnetic field lines in a region, 
we can just show the important aspects of the topology; such a picture of the magnetic 
structure is called the magnetic skeleton of the field. In 3D, we have similar properties, 
now with separatrix surfaces separating the volume into topologically distinct regions, 
and these surfaces intercept at a separator. 
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(a) (b) (c) 




Fig. 1 (a) X-type null point for = 1. This potential neutral point has separatrices (red 
lines) intersecting at an angle of tt/2. (fe) O-type null point, with a = —1. A blue star denotes 
the null point, (c) Single potential magnetic null point configuration created by interaction of 
two dipoles. Here, = y/ [(x + A)^ + y^j +y/ ^(^x - \ f + y^] , for A = 0.5. Red lines denote 
the separatrices. 



1.3.1 Mathematical description of null points 



Let us first cons ider n ull points in 2D (e.g. Dungev Il953l : [iQSSl l. Following §1.3.1 of 
Priest & Forbes (|200Cll ). we assume a magnetic field of the form: 

B^[B,,{x,y),By{x,y),0] . 

A null point occurs at the point {xq, j/q) if: 

Bx{xo,yo) = and By{xo,yo} = ■ 

Expanding Bx and By in a Taylor series about (xg, yo) gives the linear approximation; 

dBx 



Bx = 



dBx 
dx 



[x - a;o) + 

{xam) <jy 



{y - yo) 



By 



a{x - xq) + b{y ~ yo) , 
dBy , . dB, 

c(x - a;o) - a{y - yo) , 



(y - yo) 



(3) 



(4) 



where the coefficients a, b, c are arbitrary. 

Let us now introduce the vector potential (also called the fiux function). A, such 
that B = V X A, and in 2D we have A — (0, 0, Az). Thus, we have: 



B 



dAz 
dy 



dAz 
dx 



,0 



(5) 



Integrating equations Q and (|4| gives the corresponding vector potential as: 

Az = a{x - Xo){y-yo) + ^(j/- J/o)^ - ^{^ ^ xq)"^ , 

where we have chosen the arbitrary constant of integration such that Az vanishes at 
{xo,yo)- Further simplification is possible by rotating the a::y-axes through an angle 6 




Fig. 2 Magnetic fields containing two null points, (a) Magnetic configuration containing both 
X-type and O-type null points. Here = X^x - - (x - \)^/3 (B = [-2y, {x - \f - A^] ) 
for A = 0.5. Red lines/blue star denotes the separatrices/O-type null point. (6) Potential 
magnetic configuration containing two X-type null points connected by a separator. Here, 
Az = —x^y + + A^j/, where A = 1. (fe) Potential magnetic configuration containing two 
X-type null points not connected by a separator. Here, Az = —xy^ +x^ /3 — X^x, where A = 1. 
Red lines denote the separatrices. 



to give new x' , j/'-axes, and choosing the angle 6 such that tan 2^ : 
simplification gives the corresponding vector potential as: 



2L 



[y -yo 



where 



B 



2a^ + - 
^4a2 + (6 + c)2 



,2 2 

- a 



2 

a = 



Co) 



-2a/{b + c). This 



(6) 



4a^ 



2a2 + fo2 



Here, B is the characteristic strength of the magnetic field and L is the characteristic 
length-scale over which the field varies. The corresponding field components are: 



Bx 



B 



(y' 



ya) 



and B 



y 



2,1 I \ 

a. [x ~xq) 



(7) 



Magnetic field lines are defined by Az equal to a constant. For a > 0, the fieldlines 
are hyperbolic, giving an X-type null point. The separatrices are given by y' — Vq = 
iiza{x' — Xq) and are inclined at an angle iarctana to the a;'-axis. The magnetic 
fieldlines for = 1, Xq = i/q = can be seen in Figure [1^. 

The value of a (and thus the angle between the separatrices) is related to the 
current density. The current density is given by: 



j = i (V X B) = --V^AzZ = — ^ 



(l-a2)z 



Thus, a null point is potential if a = ±1 (i.e. an X-type configuration of rectangular 
hyperbola) and the angle between the separatrices is 7r/2. Note that for an O-point, 
a is imaginary and so the current density is always non-zero. Thus, O-type neutral 
points can never be potential. 

An O-type null point magnetic configuration can be seen in Figure[TjD, for = —1, 
x'q = y'fj = 0. However, note that the simple 2D magnetic field configuration of equa- 
tion is only valid close to the null point: as x' and/or y' get very large, B becomes 
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(a) (b) 




X 



Fig. 3 (a) Proper radial null point, described by B = {x,y, —2z), i.e. e = 1. (6) Improper 
radial null point, described by B = {x, ey, — [e + 1] z), for e = 1/2. Note for e = 1/2, the field 
lines rapidly curve such that they run parallel to the a;— axis along y = 0. In both figures, the 
j;— axis indicates the spine and the xy— plane at z = denotes the fan-plane. The red fieldlines 
have been tracked from the z = I plane, the blue from z = —1. 



unphysically large. Figure [TJ; denotes a more realistic single magnetic null point config- 
uration created by the interaction of two dipoles. This configuration comprises of four 
separatrices and an X-point, and as x' and/or y' get large, the field strength becomes 
smaller (i.e. a more physical field). 

Magnetic configurations can also contain multiple null points, and it can be argued 
that null points appear in pairs; a double null p oint m ay arise as a loc al bif urcation 
of a single 2D null point (see e.g. Galsgaard et a/. 1 19961 : Brown & Priest [l99j. Figure 
[2^. shows a magnetic configuration containing both X-type and 0-type null points. 
Figures [2J3 and [2J; present magnetic configuration containing two X-type null points 
connected by a separator and not connected by a separator, respectively. 



1.3.2 Three-dimensional magnetic null points 



Magnetic null points also exist in three dimensions, but occur in a different form to 
those described in i31.3.1l In 3D, potential null points are of the form: 

B = j{x,ey,-[e + l]z) , (8) 

where the parameter e is related to th e pre dominate direction of alignment of the 
fieldlines in the fan plane. Parnell et al. (|l996l ) investigated and classified the different 
types of linear magnetic null points that can exist (our e parameter is called p in their 
work). Topologically, this 3D null consists of two key parts: the 2:— axis represents a 
special, isolated fieldline called the spme which approaches the null from above and 
below (as found by Priest & Titov Il996l ) and the xy— plane through 2 = is known 
as the fan-plane and consists of a surface of fieldlines spreading out radially from the 
null. Figure |3] shows two ex amples of 3D null points: e — 1 (Figure [3^) and e = 1/2 
(Figure EJd). Titov & Hornig (|2000l ) have investigated the steady state structures of 3D 
magnetic null points. 
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Equation (|8|) is the gener al exp ression for the linear field about a potential magnetic 
nuU point (see ParneU et al. If 996l : §/!/). For e > 0, 3D nulls are described as positive 
nulls, I. e. the spine points into the null and the field lines in the fan are directed away. 
In addition, all potential nulls are designated radial, i.e. there is no spiral motions in 
the fan-plane. 

In this review paper, we only consider positive, potential null points, and thus there 
are three cases to consider: 

— e = f : describes a proper null (Figure |3^). This magnetic null has cylindrical sym- 
metry about the spine axis (so is actually only a 2.5D null point). 

— e > 0, e 7^ f : describes an improper null (Figure |3}d). Field lines rapidly curve such 
that they run parallel to the a;— axis if < e < f and parallel to the y— axis if e > f . 

— e = 0: equation ((8} reduces to a simple 2D X-point potential field in the ccz— plane 
and forms a null line along the j/— axis through x = z — 0. 



1.4 Statistics of coronal null points 

We have provided a mathematical description of null points, but how common are 
null points in the corona? Null points are an inevitable consequence of the distributed 
isolated magnetic flux sources at the photospheric surface. Using photospheric mag- 
netograms to provide the field distribution on the lower boundary, both potential and 
non-potential (nonlinear force-free) fleld extrapolations suggest that there are always 
likely to be null points in the corona. The number of such singular points will depend 
upon the magnetic complexity of the photospheric flux distribution. Detailed investi- 



gations of the coronal m agneti c field, using such po tentia l field calculations, can be 
found in Beveridge et al. (|2002l ) and Brown & Priest ( 200 ih . The properties of coronal 



null points have also been c onside red through theo retical considerations (e.g. P arnel l 
et al. 1 19961 : Brown & Priest I2OO1I : Beveridge et al. I2OO2I : Parnell & Galsgaard l2004l : 
ParneU et aZ. bOOSh . 

The statistics of coronal null points has been investigated using two metho dologies: 
direc t mea surement from potential field extrapolations (e.g. Close et al. l2004l : Regnier 
et al. I2OO8I ) and, sec ondly^ as an estimat e from the Fourier spectrum of magnetograms 
(Longcope & ParneU bOOg '). Close et al. \2004 ) calculated a potential field extrapolation 



from a high resolution MDI magnetog ram a nd found 1.7 x 10 ^ magnetic null points 



per square megameter. Regnier et al. 1 2008h performed a similar investigation using a 
magnetogram from the Narrow band Filter Imager onboard Hinode and found 6.7x 10~^ 
Mm~^. Longcope & Parnell (|2009l ) investigated 562 MDI magnetograms using the 
Fourier spectrum of magnetograms and found 3.1 x 10"'^ ± 3.0 x 10^^ coronal null 
points per square megameter (at altitudes greater than 1.5 Mm). Alternatively, we can 
estimate the total number of coronal null points by multiplying these results by the 
surface area of the Sun (i.e. to provide a crude estimate, where we assume the Sun 
is free of act ive re gions and coronal holes). This corre sponds to approximately 10, 000 
(Close et aZ. l2004h . 19,000 (Longcope & ParneU [20oi ) or 40,000 (Regnier et al. I2OO8I I 
coronal null points. 



More recently. Cook et al. (|2009l ) investigated the solar cycle variation of coronal 



null points using a potential field source surface model in spherical geometry, and find 
that there is no significant variation in the number of nulls found from the rising to 
the declining phase (indicating that null points are present throughout both phases of 
the solar cycle). 
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Several investigations al so con sider specific examples of null points in the corona. 
For example, Aulanier et al. ( 2000l 'l investigated a class M3 flare that occured on 14 July 
1998 above a 5— spot. Using potential fleld extrapolations, the authors recreated the 
pre-flare magnetic topology from Kitt Peak line-of-sight magnetograms and reve alled a 
single coronal null point located above the 5— spot. Secondly, Ugarte-Urra et al. { 200?! ') 
investigated the magnetic topology of 26 CME events by performing potential field 
extrapolations from corresponding MDI magnetograms, and find that magnetic null 
points are present in a large number of the pre-CME topologies. 

Finally, we note t hat a null point plays a key role in the magn etic breakout mode l 
(e.g. Antio chos I l998l : Antiochos et a/. 1 19991 : MacNeice et al. l2004l : Lynch et al. l2004l : 



Choe et The equilibrium set-up of the magnetic breakout model model con- 

sists of a quadrupolar photospheric fiux distribution coupled with an overlying field, 
and such a set-up contains a coronal null point. Such a null point is an ideal candidate 
for the study of MHD wave behaviour about coronal null points (i.e. the investigations 
detailed in this review paper). However, it is important to stress that the null point 
in the breakout model is not the only candidate - the ideas and investigations detailed 
below apply equally well around coronal null points found elsewhere (i.e. quiet Sun and 
inside active regions). Thus, the coronal null points we are describing in this paper are 
not solely those involved in the magnetic breakout model. 



1.5 Why is this area of study interesting or important? 

The motivation for investigating the behaviour of MHD wave propagation in the neigh- 
bourhood of magnetic null points can be summarised as follows: 

— MHD wave propagation in inhomogeneous media is a fundamental plasma process, 
and the study of MHD wave behaviour in the neighbourhood of magnetic null 
points directly contributes to this area. 

— We now know that MHD wave perturbations are omnipresent in the corona. We 
also know that null points are an inevitable consequence of the distributed isolated 
magnetic flux sources at the photospheric surface (where the number of such sin- 
gular points will depend upon the magnetic complexity of the photospheric flux 
distribution). Thus, these two areas of scientific study; MHD waves and magnetic 
topology, will encounter each other at some point, i.e. MHD waves will propagate 
into the neighbourhood of coronal null points (e.g. blast waves from a flare will at 
some point encounter a null point). Thus, the study of MHD waves around null 
points is itself a fundamental coronal process. 

— The study of MHD wave behaviour in the neighbourhood of magnetic null points is 
also interesting in its own right and, as we shall see, often provides critical insights 
into other areas of plasma physics, including: mode-conversion f ijS.SII . reconnection 
(!j3]6l), quasi-periodic pulsations f i)3.7|) and phase-mixing (53]8|. 



2 Two-dimensional null points in a cylindrically symmetric geometry 

The flrst investigation into the behaviour of MHD waves in the neighb ourhood of 2D 
magnetic null point was performed by Bulanov & Syrovatskii (198Cj). The authors 
considered the linearised MHD equations (equations [T| under the cold plasma approx- 
imation (/3 = 0) in a cylindrically symmetric geometry, where a circular boundary is 
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imposed at r = 1 at which the fieldlines are held fixed. Their investigation produces a 
detailed discussion of the propagation of fast and Alfven waves in the neigbourhood of 
an equilibrium magnetic field: 

Bo = f (-y,-:r) , (9) 

which corresponds to the magnetic configuration in Figure [1^. The subscript denotes 
an equilibrium quantity. 

Bulanov & Syrovatskii also noted that in this 2D geometry, the z— component mo- 
tions are decoupled from xj/— plane motions. Thus, the Alfven wave, which is governed 
by motions transverse to the magnetic field (i.e. Vz z) is decoupled from the magnetoa- 
coustic waves (motions in the sy— plane). Hence, it is possible to consider the Alfven 
wave and the magnetoacoustic waves separately in a 2D geometry. 

In their paper, harmonic fast waves are generated at the r = 1 boundary and these 
propagate inwards towards the null point, and Bulanov & Syrovatskii find that, in the 
asymptotic limit r — >■ 0, the perturbations have azimuthal symmetry, i.e. propagate 
as cylindrical waves. This was the first piece of work that indicated a key relationship 
between fast waves and null points. However, the assumed cylindrical symmetry means 
that the disturbances can only propagate either towards or away from the null point, 
and are already encircling the null point. Thus, it is unclear if this is a general result. 

To investigate the Alfven wave, Bulanov & Syrovatskii make a coordinate transfor- 
mation such that the coordinate lines coincide with the lines of force, where: 

C = i (a;^ - y^) = ^r^ cos 261 ,rj = xy = ^r^ sin 261 . (10) 

The authors find that Alfven perturbations propagate along magnetic fieldlines at 
the local Alfven speed. The inhomogeneity of the Alfven speed profile leads to an 
exponential increase in the gradients in the system, and in the asymptotic limit t — > oo, 
these gradients accumulate at the separatrices of the field, i.e. ( — 0. Again, this 
result was the first to indicate a relationship between Alfven wave propagation and the 
location of the separatrices. 

An alternative approach to the study of MHD wave behaviour in the neighbourhood 
of a 2D magnetic null point was investig ated in a series of pap ers by Craig and co- 
workers (Craig & McClvmont ll99ll : 11993 : Craig & Watson [l99a ) in which the authors 
considered perturbations of the flux function Az- Here, the focus was on applications for 
magnetic reconnection, specifically to investigate if null points (viewed as weaknesses 
in the magnetic field) could collapse in response to imposed boundary motions. To 
clearly demonstrate their results, we repeat part of their analysis here: 

In polar coordinates, an X-type null point, located at the origin, can be expressed 

as: 

Bo = — (r sin 26) f -f r cos 261 ^) , (11) 
it 

where Bq denotes the equilibrium magnetic field and R represents the typical size of 
a coronal magnetic structure. 

The MHD equations ^ can be simplified by invoking the flux function Az (equation 
[S]). In polar coordinates, the equilibrium flux function (^lo) for ^ simple 2D X-point 
(i.e equation (SI is: 



An = cos 261 

^ 2 



(12) 
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where the evolution of Az is governed by the induction equation and the dynamics are 
governed by the momentum equation: 

^ + (v ■ V) A, = 77V"A. , 
at 

These equations can be linearised about the equilibrium flux function such that Az = 
Aq + Ai, and combined to form a single differential equation for the perturbed flux 
function (Ai): 



i = lVAopv2^i+77V2^ , (13) 



dt^ ' ' ' dt 

where, from equation H12II . |VAoP = ■ 

We can see that the right-hand-side of equation (|13|l has two parts. Consider a 
region close to the origin, Tc, where Tc ~ rj^^"^ is the skin depth. Far from the origin, 
r ^ Tc, advection effects dominate and equation (|13|) reduces to a wave equation: 



d^A 



at2 



i = r^V^Ai (14) 



Here, the rate of propagation of information is governed by the wave speed proportional 
to r which makes the signal travel time logarithmic in r. Thus, a disturbance on the 
outer boundary (r = 1) propagates into the diffusion region in a time 5t ~ | Inrd ~ 
\\nr). 

Alternatively, close to the origin, r ^ rc, and equation (|13|l reduces to the diffusion 
equation: 

Diffusion is ultimately responsible for dissipating the kinetic and magnetic energy in 
the system. 

Next, let us assume separation of variables such that: 
Ai(r,e,t) = /(r)e™V* , 



where m is the azimuthal wavenumber. The eigenequation for /(r) is then: 

(15) 



dr \ dr J 



(l + ryA/r2) 



- m 



Close to the origin (r <^ Tq), the radial dependence of the current jz ~ —V'^Ai = 



J. 



(-AA?) 



1/2 



We note that only the Jg Bessel function is nonzero at the origin 

(which is equivalent to a nonvanishing displacement in the flux function) and so the 
m = mode is the only mode which allows topological reconnection. Assuming m = Q, 
the evolution of Az can then be obtained from the analytical or numerical solution of 
equation (|13|) . 



Craig & McClymont (|l99ir i derived equation (|13p and investigated the fieldlines 
passing through the (circular) boundary in a particular manner in order to perturb the 
field - shifting the footpoints so as to "close up" the angle of the X-point. The resulting 
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field perturbations cause tiie null point to collapse to form a current sheet in which 
reconnection can release magnetic energy. In these models the boundary motions move 
the field lines but do not return them to their original positions (akin to photospheric 
footpoint motion). Thus, the Poynting flux induced by the imposed motion (and then 
fixing the fleld after the motion is complete) accumulates at the resulting current sheet 
and provides the energy released in the reconnection. 

Craig & McClymont find that, as their system evolves, the fieldlines reconnect as 
they pass through the null point (located at the origin), indicating that resistivity is 
essential to this mode. Their initial vertical 'current sheet' begins reconnecting and 
the inertia of the flowing plasma carries the system past the equilibrium configuration, 
until a weaker horizontal current sheet is formed. Then, a much weaker vertical current 
sheet returns at a later time (one complete cycle). After three such cycles, the system 
is close to its equilibrium (potential) configuration. The reconnection is found to be 
oscillatory, with inertial overshoot of the plasma carrying more magnetic flux through 
the neutral point than is required to reach a static equilibrium. The reconnection rate 
scales as |ln77p, and so the reconnection is described as fast. This is an important 
results as it shows that the damping of the fast wave is still highly significant even 
when the resistivity is small. 

Craig & Watson (|l992l ) considered the radial propagation of the m — Q mode 



and solved equation H15|) using a mixture of analytical and numerical solutions. They 
demonstrated that the propagation of the m = wave towards the null point generates 
an exponentially large increase in the current density and that magnetic resistivity dis- 
sipates this current in a time related to log 77, in agreement with Craig & McClymont 
(|l99lh . Their initial disturbance is given as a function of radius, i.e. an internal per- 
turbation is considered. In their investigation, the outer radial boundary is held fixed 
so that any outgoing waves will be refiected back towards the null point. This means 
that all the energy in th e wave motions is contained within a fixed region. 

Craig & McClymont (|l993l ) investigated the normal mode solutions for both m = Q 
and m 7^ modes with resistivity included. Again they emphasised that the current 
builds up as the inverse square of the radial distance from the null point. Craig & 
McClymont ( 19931 ) also explicitly report on the focusing of the wave (energy) onto the 



neutral point due to the gradient in Alfven speed (clearly seen in equation [14]) . 

Independently, Hassam ((1992 ') also investigated the behaviour of stressed X-points 
in a cylindrically symmetric geometry. Hassam performed a similar derivation to that 
of equation (|13p and recogni sed th at equation (|15p can be recast as a hypergeometric 
equation (e.g. Oberhettinger ^9901 ) when m = 0: 

d / d A /A^ ^ 



using the transformation z = —r'^/rjX. The relaxation time found using thi s formulatio n 
is I Inryl^ and so i s in a greement with that found by Craig & McClymont (|l99ll : ll993l '). 



Craig et al. (|2005l l comment that it is surprising that fluid viscosity has been 
neglected in previous studies, and state that the leading terms in the viscous stress 
tensor actually dominate the plasma resistiv ity by many orders of magnitude for typical 
coronal plasmas (as e mphas ised by Hollweg |l986h . These authors extend the model of 



Craig & McClymont ( 19931 1 to include the effect of (isotropic) scalar fluid viscosity. 



They find the inclusion of viscosity can have a dramatic effect: for non-reconnective 
modes (m > 0), the dissipation rate depends only logarithmically on the magnitude 
of the dissipative coefficient (thus, for rj > f the problem is dominated by ln?7, and 
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> 77 it is dominated by \nv). For reconnective disturbances (?n — 0), the oscillatory 
reconnection is suppressed if v > rj, whereas for rj > v, fast oscillatory reconnection 
is regained. However, the authors do note that although viscosity can dramatically 
influence the rate of dissipati on in the system, finite resistivity is still required for 
reconnection to occ ur. Cr aig l|2008l ') has extended this study to include non- isotropic 
viscosity (Braeinskii I igGSl ) and find that the main results of Craig et al. (|2005h are still 
valid. 

These papers have led the way in understanding MHD motions in the neighbour- 
hood of a 2D X-point. However, the assumed cylindrical symmetry means that the 
magnetoacoustic disturbances can only propagate either towards or away from the null 
point, and attention has been restricted to a circular reflecting boundary, so all outgo- 
ing waves are reflected back into the vicinity of the null point. Thus, in a sense there 
is nowhere else for the wave to propagate except into the null poi nt. 

In addition, in all these papers (except Bulanov & Svrovatskii 1 198ol which consid- 
ered asymptotic limits) the boundary motions move the field lines but do not return 
them to their original positions. The Poynting flux induced by this imposed motion 
provides the energy released in the resultant current sheet. However, if the boundary 
motions are simply due to the passing of incoming waves through the boundary, then 
it is not clear that the null point need collapse and form a current sheet. Furthermore, 
if this is the case, then it is not clear if the energy in the wave (again due to the 
Poynting flux through the boundary) will dissipate or simply pass through one of the 
other boundaries. 

Finally, all these papers have assumed that the system is best described in terms 
of normal modes, where a single normal mode can be thought of as the long time 
evolution of a system. However, normal mode analysis does not allow us to concentrate 
on the transient features of the wave propagation. 

For these reasons, it is informative to specifically track the propagation of boundary- 
driven (as opposed to internally generated) asymmetric disturbances into the domain. 



3 Two-dimensional null points in a cartesian geometry 

The first investigation that specifi cally made the exten sion a way from polar symmetry 
was performed by Ofman j 19921 ) and Ofman et al. (|l993h . who performed nonlin- 



ear, resistive 2D MHD calculations of the reconnection in the stressed null point and 
obtained the |lnr;|^ scaling law numerically, as well as solving the linear dispersion 
relation (equation [TS]) for all azimuthally nonsymmetric pert urbations (m > 0) analyt- 
ically in agreement with the results of Craig & McClymont (|l99ll : ll993l ) and Hassam 



(|1992| ). However, the authors also stressed the significance of the choice of boundary 
condit ion, a nd the strong influence they have on the permissible solutions. Steinolfson 
et al. (|l995l ') investigated the effects of the boundary conditions further by perfoming 
nonlinear, resistive MHD simulations to study stressed X-points, contrasting both rigid 
and open boundary conditions. The authors found tha t, for r igid b oundary conditions. 



they could recover the results of Craig & McClymont (|l99ll : ll993l l. but found that for 



open boundary conditions, the X-point was instead deformed by their perturbation to 
form a current sheet. 



Shortly thereafter, Hassam & Lambert (|l996l ) utilised a cartesian geometry to in- 
vestigate the propagation of the Alfven wave in the neighbourhood of a simple X-point 
(of the form seen in Figure [1^). Using the same coordinate transformations (equation 
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I10|l as in Bulanov & Syrovatskii ([1980|), Hassam & Lambert again find that the Alfven 
wave propagates along magnetic fieldhnes, i.e. the fluid elements are conflned to the 
magnetic fieldlines they are generated on. However, Hassam & Lambert also stress 
the distinction between two different types of boundary driven Alfven wave motions. 
Perturbations that initially straddle the separatrices eventually accumulate along the 
separatrices, whereas disturbances that do not initially straddle the separatrices do not 
accumulate along the separatrices, but l eave the system and are eventually damped by 
phase-mixing (Heyvaerts & Priest 1 1983h . 

Hassam & Lambert investigated the propagation of the Alfven wave in a square 
numerical box by driving a harmonic wave train (polarised transverse to the plane of 
the magnetic field) on both the left and right boundaries, simulating periodic footpoint 
motion. Their choice to drive wave-trains into the domain from opposite sides of the 
box means the forcing is antisymmetric in the horizontal direction. In addition, line- 
tied conditions are used on the upper and lower boundaries, i.e. both velocity and the 
normal derivative of magnetic field are kept zero. Thus, wave motions are reflected 
back into the numerical domain, limiting the evolution. 

An investigation of MHD wave propagation in a /3 = plasma in the neighbourhood 
of a 2D nul l point has been looked at in a series of papers by McLaughlin & Hood 
( 20041 : 12005l : l2006al ). where the focus was on more general disturbances, more general 
boundary conditions and single wave pulses (rather than harmonic wave trains). By 
looking at boundary-driven disturbances generated from a single boundary, coupled 
with non-reflecting boundary conditions, such investigation allow us to focus on the 
transient features on the propagation. 

To clearly demonstrate the results of McLaughlin & Hood, we repeat part of their 
analysis here: 



3.1 Governing equations and coordinate system of McLaughlin & Hood 

To study the nature of wave propagation near null points, McLaughlin & Hood utilised 
the linearised MHD equations. Here, we use subscripts of for equilibrium quantities 
and 1 for perturbed quantities, such that v = -|- vi (no background flows). The only 
exception is B = Bq + b, where b = {bx, by,bz). 

McLaughlin & Hood now consider a special coordinate system for vi: 











|Boi 




|Bo| 





where Aq is the equilibrium vector potential. The terms in brackets are unit vectors. 
This splits the velocity into parallel and perpendicular components. This will make our 
MHD mode detection and interpretation easier. For example, in a low /3— plasma, the 
slow wave is guided by the magnetic field and has a velocity component that is mainly 
field-aligned. This makes perfect sense when /3 ^ 1 but its usefulness is less clear when 
/? ^ 1 near the null point (recall the properties detailed in Table 1). 

Now consider a change of scale to non-dimensionalise; let vi — vv^, vj^ = vBy*^, 
V|| = vBvJJ, Bo = -BBq, b = Bh*, x = Lx* , z = Lz* , pi = poPi, V = V*/L, t = it*, 
Aq — BLAq and r] = rjo, where we let * denote a dimensionless quantity and v, B, L, 
Pq, t and 770 are constants with the dimensions of the variable they are scaling. We then 
set B I JJTpQ — V and y = L/t. We also set rj^t/Lp' = , where Rm is the magnetic 
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Reynolds number, and set (Sq = 2iipo/B^, where /?o is the plasma-/? at a distance unity 
from the null/origin. 

This process non-dimensionalises the hnearised MHD equations. For the rest of this 
section, we drop the star indices; the fact that the variables are now non-dimensionahsed 
is understood. Thus, the linearised, non-dimensionahsed equations are: 

Po^vi = -IBol' [(V X b) ■ y] + ^VAo ■ Vpi , 



Po% = (Bo-V)&. 



dt 
dhx 

dt 
dby 

dt 

dt 

dpi 
dt 



V > 

[(Vv_L X y) • + -^V^bx , 

rim 

(Bo- V)«y + -^v2fe„ , 

V hz , 



= [(Vv^ X y) ■ z] 



Bovii 



R. 



V I \/A, 



(16) 



Note that the 2D geometry considered by M cLaughlin & Hood is in the xz— plane. 
As noted first by Bulanov & Syrovatskii l|l980l ). this means that the Alfven wave (in 
the y— direction) is decoupled from the magnetoacoustic waves (in the 2:2;— plane). 

Finally, we note that equations p6|) simplify greatly if we neglect pressure pertur- 
bations (i.e. assume /3 = 0). Under the /3 = assumption, the plasma pressure plays no 
part in the dynamics of the system, and so the linearised equation of mass continuity 
has no influence on the momen tum e quation and so in effect the plasma is arbitrarily 
compressible (Craig & Watson Il992l ) and we assume the background gas density is 
unifo rm (po). A spatial varia tion in po c an ca use phase mixing (Heyvaerts & Priest 
Il983l : De Moortel et a/. 1 19991 : Hood et al. hooi ). 

The ideal, P — 0, linearised MHD equations naturally decouple into two equations 
for the fast magnetoacoustic wave (governed here by vj^) and for the Alfven wave 
(governed by Vy). The slow wave is absent in the 13 — limit (vy =0). 

Under these assumptions, the linearised equations for the fast magnetoacoustic 
wave can be combined to form a single wave equation: 



d^ 
dt'^ 



2 

v_l 



+ 



d' 



\ dx^ d. 



(17) 



where va(x,z) = |Bo|/y^po = y [Bi + Bij /pQ is the equilibrium (unperturbed) 
Alfven speed. 

Similarly, the linearised equations for the Alfven wave can be combined to form a 
single wave equation: 



§^^y-iBo-Vfvy = {Bx^^+Bz§- 



(18) 



Wave equations H17|) and H18|) are the primary equations governing the behaviour 
of the linear, P — 0, fast and Alfven waves in an equilibrium magnetic field Bg. These 
equa t ions fo rm the basis of the investigations carried out by McLaughlin & Hood (|2004l : 
I2005l : l2006al ) for various magnetic configurations. 
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Fig. 4 (a) Contour of at t = 2 for a fast wave driven from the upper boundary. The 
overplotted black lines are the WKB solution, where the three lines represent the leading, 
middle and trailing edges of the wave pulse. The blue cross denotes the null point (located 
at the origin), (b) Ray paths of the WKB solution for an Alfven wave driven at the upper 
boundary, after a time t = 2, for starting points of a; = 0, 0.025, ...0.5. 



3.2 Single two-dimensional null point 



McLaughlin & Hood (|2004l ) investigated the behaviour of the fast and Alfven waves 
about a simple 2D X-type null point using the following equilibrium magnetic field: 



Bo = 



(19) 



This magnetic field represents a 7r/4 rotation of the magnetic field seen in Figure [TJi. 
Note that this particular choice of magnetic field is only valid in the neighbourhood of 
the null point located at g = , 2 = 0. 

McLaughlin & Hood l |-'IHI|| l considered a single wave pulse coming in from the top 
boundary of the form: 



Vj_(x-,2max) = 



sin Ljt for Q<t< ^ 
otherwise 



9vx 
dx 



= 



dx 



= 



(20) 



dz 



= 



where uj is the frequency. The authors find that the linear fast magnetoacoustic wave 
travels towards the neighbourhood of the X-point and bends around it. Since the Alfven 
speed, Vji{x,z) = jBol/y^po = x'^ + z'^ , is spatially varying, the wave travels faster the 



further it is away from the null point. Thus, the wave demonstrates refraction and this 
can be seen in F igure 2^. A similar refraction phenomenon was found by Nakariakov 
& Roberts (|l995l ). This refraction effect wraps the wave around the null point, and it 
is this that is the key feature of linear, /3 = fast wave propagation. This refraction 
effect is the (non-radial) g eneral isation of the (purely radia l) focusing effect reported 
by Bulanov & Syrovatskii (|l980l ) and Craig & McClymont (|l99ll : Il993l ). 
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Note, since the Alfven speed drops to zero at the null point, the wave never reaches 
there, but the length scales (which can be thought of as the distance between the leading 
and trailing edges of the wave pulse) rapidly decrease, indicating th at gra dients, and 
hence the current density, will rapidly increase. McLaughlin & Ifood (|2004l ) show that 
for the simple 2D X-point, all gradients increase exponentially as they approach the 
null point. The rate of this build-up is extremely important, as it implies that resistive 
dissipation will eventually become important, regardless of the size of rj, and this will 
convert the wave energy into (ohmic) heat. In fact, the exponential growth of the 
current density indicates that the time for magnetic diffusion to become important 
will depend on In 77. 

This is in good agreement with Craig & McClymont (|l99ll : I1993I ') and Craig & 
Watson (|l992l ) who had previously found that the reconnection rate scales as | lnry|^. 
This means that wave dissipation will be very efficient, and predicts that null points 
will be the natural locations of linear fast wave energy deposition and preferential 
heating. 

To confirm their results, McLaughlin & Hood ( 2004[ ) also solved eq uation ((T7| 
approximately using the WKB approximation (e.g. Murrav fl927l : Sneddon [l957l '). The 
WKB solution of McLaughlin & Hood as sumes u ^ 1 and is obtained using the method 
of characteristics (e.g. Bender & Orszag ri978h and is in excellent agreement with the 
original numerical solution. This can be seen from the overplot in Figure 3^. 

By considering equation HI8II . McLaughlin & Hood (|2004l ^ find that the linear 
Alfven wave propagates down from the upper boundary and be gins t o spread out, 
following the field lines, in agreement with Bulanov & Syrovatskii ( 19801 ). The wave is 
confined to the fieldlines it is excited on. As the wave approaches the separatrix (defined 
by the a::— ax;is), the pulse thins but keeps its original amplitude. The wave eventually 
accumulates along the separatrices. As for the fast wave, we have decreasing length 
scales, and for this choice of set-up, jx grows exponentially in time. Hence, the authors 
find that the Alfven wave causes curre nt den sity to accumulate along the separatrices 
(in agreement with Hassam & Lambert ll996l l. Hence, all the Alfven wave energy will be 
dissipated along the separatrices and these will be the locations for preferential heating. 

A WKB solution was also obtained for the Alfven wave and this was in excellent 
agreement with the numerical results. In Figure [JIj, we can see the evolution of fluid 
elements that begin at points x = 0,0.025, ...0.5, which clearly demonstrates how the 
fluid elements simply travel along the fieldlines they start on. In addition, note that 
at t = 2 the fluid elements have all travelled different distances along their respective 
fieldlines, but the wave remains planar. McLaughlin & Hood explain this and show 
that the leading edge has in fact reached z = Znia.xe "* = 2e"2 = 0.27. 

Note that the behaviour of the Alfven wave is different to that of the fast wave in 
the sense that the two wave types deposit all their wave energy at different areas (along 
separatrices as opposed to at the null point). However, the phenomenon of depositing 
wave energy in a specific area is common to both. In addition, this work is in good 
agreement with previous work in cylindrically symmetric geometries {^^. Thus, some 
general and robust properties of the fast and Alfven wave are becoming apparent. 



3.3 Pair of two-dimensional null points 



McLaughlin & Hood (|2005l ) repeated the aboce investigation for a magnetic field con- 
taining two null points. Such a configuration is particularly relevant since it can be 
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argued that nuU points appear in pairs. For example, a double null po int m ay arise as 
a loca l bifu rcation of a single 2D null point (see e.g. Galsgaard et al. Il996l : Brown & 
Priest Il998h . 

Two magnetic configurations were considered; one containing a separator and one 
that does not. The magnetic configuration with a separator is taken as: 

Bo^^{x^-z^-\^,0,-2xz) , (21) 

where 2A is the distance between the null points. Note that this introduces a char- 
acteristic length scale into the system, whereas previously the single X-point had no 
charactistic length scale. The other equilibrium magnetic field considered takes the 
form: 

Bo = (2x2,0,2:^ -2^ -A^) . (22) 

These equilibrium magnetic field configurations can be seen in Figures [2b andl^J:. The 
authors considered the /3 = linearised MHD equations and solved equations H17|) and 
(|18|1 for these two magnetic equilibria. 

McLaughlin & Hood (|2005h find that for the linear fast magnetoacoustic wave 



approaching the two null points from above, the wave travels down towards the null 
points and begins to refract around them both. The wave the 'breaks' into two along 
the line x = (due to symmetry), with each half of the wave going to its closest null 
point. Each part of the wave then continues to wrap around its respective null point 
repeatedly, eventually accumulating at that specific null point. 

In the case of the fast wave pulse travelling in from the side boundary, we see a 
similar effect (i.e. a refraction effect, wave breakage and accumulation at the nulls), but 
in this case the wave is not equally shared between the null points. For example, for a 
fast wave travelling in from the left boundary, initially the pulse thins, begins to feel 
the effect of the left-hand-side null point and begins to refract around this null. As the 
ends of the wave wrap around behind the left null point, they then become influenced 
by the right-hand-side null point. These arms of the wave then proceed to wrap around 
the right null point, flattening the wave. Furthermore, the two parts of the wave now 
travelling through the area between the null points have non-zero Alfven speed, and so 
can propagate through this area. These parts of the wave break along a: = and then 
proceed to wrap around the null point closest to them. 

As before, it is clear the refraction effect focu ses all the energy of the incident wave 
towards the null points, but McLaughlin & Hood (|2005h find that the angle that the fast 
wave approaches the null points from will determine what proportion of wave energy 
ends up at each null point (i.e. where the wave 'breaks' ). In the case of the fast wave, 
all the wave energy is accumulating at the null points and since we have a changing 
perturbed magnetic field with increasing gradients, this is where current density will 
accumulate. Thus, fast wave heating will naturally occur at both null points. 

For the Alfven wave, the results show that the wave propagates along the field 
lines, thins but keeps its original amplitude, and eventually accumulates along the 
separatrices (again, in agreement with previous studies). 

Th us, M cLaughlin & Hood (|2005l ') find that the key results from McLaughlin & 
Hood (|2004l ) carry over from a single 2D null configuration to that of a pair of 2D null 
points. 
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3.4 Single null point configuration created by two magnetic dipoles 

McLaughlin & Hood (|2006al ') again investigate the behaviour of the fast and Alfven 
wave in an ideal, /3 = plasma, but now address a key problem with the first two 
papers: that the simple null points considered (i.e. equations 1191 [211 122|) are only valid 
locally, because as x and z get very large, Bq als o gets unphysically large. 



To address this issue, McLaughlin & Hood (|2006al 'l investigate the behaviour of 



MHD waves near an equilibrium magnetic field created by two dipoles: 

(x + Xf+z"^ -{x-Xf + z'^\ 



[{x + Xf + z-^] [[x-Xf+z-^] 



\[{x + xf + z^f [(x-A)V^2]V 

where 2A is the separation of the dipoles. This magnetic field (with 2A = 1) can be seen 
in Figure[l};. It comprises of four separatrices and an X-point located at {x, z) = (0, A). 
Note that as a: or z gets very large, the field strength becomes small. Hence, this is a 
more physical field than those previously investigated. 



McLaughlin & Hood (l2006a|) drive a fast-wave planar pulse on their lower boundary 



(along z = 0). As the wave propagates upwards, the planar wave is distorted by the 
two regions of high Alfven speed and the wave forms two peaks with maxima located 
over the loci of the magnetic field. Close to the null point, the fast wave begins to 
refract around the null point. Meanwhile, the rest of the wave (referred to as the 
wings) continue to propagate upwards and spread out (since the fast wave propagates 
roughly isotropically) . The wave is stretched between its two goals (part wrapping 
around the null and part travelling away from the magnetic skeleton) and this leads 
to the wave splitting; near the regions of high Alfven speed the localised high speed 
thins the wave and forces the split. Thus, part of the (now split) wave spirals into the 
null and the other part propagates away from the magnetic skeleton. Finally, a WKB 
solution demonstrates that there is a critical radius of influence within which a fast 
wave will be captured by the null point and that, for the parameters considered by 
McLaughlin & Hood, 40% of the wave is trapped by the null. This makes intuitive 
sense: if the fast wave is t oo far away from the magnetic null, it will not feel its effect. 
McLaughlin & Hood (|2006al ') also looked at the behaviour of the Alfven wave and 



found that, as before, the propagation follows the magnetic fieldlines but that now only 
part of the wave accumulates along the separatrices and the other part of the wave 
appears to propagate away from the magnetic skeleton. However, the Alfven wave is 
actually just following the fie ldline s (and these are spreading out) recovering the key 
results of Hassam & Lambert ( 19961 ). i.e. there is a key difference between the evolution 



of perturbations that initially straddle the separatrices and those that do not. 



3.5 MHD wave behaviour at /3 7^ null points 

The key results from McLaughlin & Hood (|2004l : I2OO5I : booeal l demonstrate that the 
behaviour of the fast wave in a /3 = plasma is entirely dominated by the Alfven-speed 
profile, and since the magnetic field drops to zero at the X-point, the wave will never 
reach the actual null. The next step is to extend the model to include plasma pressure 
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(/? 7^ plasma). The most obvious efltect is the introduction of slow magnetoacousUc 
waves. The fast wave can now also pass through the null point (there is a non-zero 
sound speed there) and thus perhaps take energy away from that area. Such a model 
could also involve mode coupling in areas where the sound speed and Alfven speed 
become comparable in magnitude. 

Such an investigation has be en carri ed out into the behaviour of magnetoacoustic 
waves by McLaughlin & Hood (|2006bl l for a single 2D null point (i.e. equilibrium 



magnetic field given by equation I19p . Note that the plasma pressure plays no role in 
the pr opaga tion of the linear Alfven wave, and so the description by McLaughlin & 
Hood (|2004h remains valid. 

In most parts of the corona, the plasma-/? (equation [2| is much less than unity and 
hence the pressure gradients in the plasma can be neglected. However, near null points 
the magnetic field vanishes and so the plasma-/? can become very large. Thus, under- 
standing the changing plasma-/? is of key importance here. Considering equilibrium 
quantities; 

^ ^ _2i£PqI^_ ^ ^ ^ /?o ^ /?o 

[x^ -f 2^) 2;2 -|_ 2;2 ^2 ' 

where r = x + z and /?o ~ 'ijJ.poL / B . 

Thus, the plasma-/? varies through the whole region, since magnetic field is varying 
everywhere throughout our model. In fact, the plasma—/? is infinite at the null point. 
In particular, we note that outside a radius of unity we have a low-/? environment and 
inside we have a high—/? environment. This will have important consequences as fast 
and slow waves have differing properties depending upon their environment (see Table 
1). 

There is also coupling between the perpendicular and parallel velocity components 
(when /? 7^ 0) and this coupling is most effective wher e the sound speed and the Alfven 
velocity are comparable in magnitude. Bogdan et al. ( 20031 ) call this zone the magnetic 



canopy or the /? ~ 1 layer. Note that here, we use the terminology /? for the (true, 
varying) plasma-/? and /?o for the (constant) value of the plasma-/? at a radius of unity. 

The /? = 1 layer occurs at radius r = However, it is not the /? = 1 layer 

that is most important in understanding this system, but the layer where the sound 
speed is equal to the Alfven speed, i.e. Cs = v^. Recalling that Cs = •\/7/?o/2fA! this 
means that the Cs = layer (or alternatively the /? = 2/7 layer) occurs at a radius 
r = •^7/?o/2. Of course, the difference between the /? = 1 layer at r = y^A) the 
Cs = VA layer at r = •\/ 7/?o /2 is very small, and hence it is easier to refer to the /? « 1 
layer. 

Finally, we note that the basic fast wave speed in this /? 7^ system is: 



cLt = \(yA + (?s) + \J {v\ + clf -^v\cl cos2 e 



2 

which for perpendicular propagation reduces to 



Cfast = wi + = |/?0 + -f 2^ . (24) 



McLaughlin & Hood (|2006bh solve the linearised MHD eq uation s (fT6)l with /Jq / 0. 



Using identical boundary conditions to McLaughlin & Hood (|2004h . i.e. eauation l20l a 
wave pulse is driven in the perpendicular velocity component on the upper boundary, 
which corresponds to driving a low-/? fast wave. It is found that the low-/? fast wave 
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propagates into the neighbourhood of the null point and begins to refract around it. 
However, as the wave crosses the Cs = layer, the low-/3 fast wave transforms into 
a high-/? fast wave (due to the change in environment) and also generates a high-/3 
slow wave. The fraction of the incident wave converted into slow wave is found to be 
proportional to Pq. The magnetoacoustic propagation now proceeds in three ways: 

— Firstly, the generated slow wave spreads out along the fieldlines, eventually accu- 
mulating along the separatrices. 

— Secondly, the remaining part of the fast wave inside the Cs = «a layer continues to 
refract and some of it now passes across the null. We identify this wave as a high-/? 
fast wave. The high-/? fast wave can pass through the null point because, although 
Vj^{0,0) = 0, there is now a non-zero sound speed there (clearly seen in equation 
I24|l . After it has crossed the null, the high-/3 fast wave continues to propagate 
downwards and crosses the Cs = va layer for a second time. As it emerges, the 
wave now becomes a low-/J fast wave and spreads out isotropically. 

— Finally, the (low-/3) fast wave located away from the null and away from the Cs = Vj^ 
layer (again referred to as the 'wings' of the low-/3 wave) are not affected by the 
non-zero sound speed (as Vj^^ Cs) and so here the refraction effect dominates. In 
fact, as these wings wrap around below the null point, they encounter the high-/? 
fast wave as it is emerging from the Cs = layer. This results in a complicated 
interference pattern, but it appears that the two waves pass through each other 
(due to the linear nature of the system). 

It is clear that there are two competing phenomena: (a) the refraction effect due to 
the varying Alfven speed and (6) a non-zero sound speed at the null which allows the 
fast wave to pass through it. It is the value of Pq that dictates which effect dominates. 
Thus, two extremes can occur. The first occurs when Pq — 0, in which c ase th e 
refraction effect dominates and we recover the results of McLaughlin & Hood (|2004l l. 
The second occurs when /3o — >■ oo and the system becomes hydrodynamic. In this case, 
the fast wave reduces to an acoustic wave and so completely passes through the null 
(effectively, it does not even see the magnetic field, since Vy^ <ti Cs). Thus, it is possible 
to understand the whole s pectrum of values of the parameter /?o ■ 

McLaughlin & Hood (|2006bl ) also noted a tendency for their system to develop 
several lobe-like structures in various variables. This occurs because the choice of mag- 
netic equilibrium naturally leads to a sin 20 and sin 49 dependence in v|| when vj^ is 
driven (this was demonstrated in polar coordinates, using equilibrium magnetic field 
lllfl . Thus, with this choice of null point, if we drive any of the velocity variables then 
the system will naturally develop a 6— dependence. 

Finally, the authors compared their results with an analytical WKB approximation. 
This resulted in two separate wave descriptions, corresponding to the fast or slow wave. 
However, the WKB solution could only reproduce the propagation of the fast wave 
(in both its low and high-/3 environments) or the slow wave but not both together. 
Instead, the WKB solution broke down at the conversion layer (at Cs = u^) where the 
approximation becomes degenerate. Thus, the WKB solution (in the form presented by 
McLaughlin & Hood) cannot be used to investigate mode conversion. Instead, addition 
terms are needed in the approximation (the authors only consider the first-order terms 
in the WKB approximation). Alternatively, this degen eracy can be overcome by using 
the method developed by Cairns & Lashmore-Davies (| 19831 ) to match WKB solutions 



across the mode conversion layer. This has been done in ID (McDougall & Hood l2007l ) 
but a 2D investigation has yet to be completed. 
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3.6 Nonlinear simulations in the neighbourhood of a 2D X-type null point 



McLaughlin & Hood (|2004l : l2005l : l2006al : l2006bh investigated the behaviour of the 



linear fast and slow magnetoacoustic waves and Alfven waves in the neighbourhood of a 
variety of 2D null points using equations (|16|l . However, the validity of the linearisation 
is questionable once the perturbed velocit y beco mes comparable to the magnitude of 
the loc al Alf ven speed. McLaughlin et al. ([20091 ') extend the models of McLaughlin & 
Hood (|2004h and Craig & McClvmont lligQlh to include nonlinear effects in a /? / 
plasma, and consider the behaviour of the nonlinear fast wave. The authors solved the 
nonlinear, compressible, resistive MHD equations (e quatio ns [T]) using a Lagrangian- 
remap, shock-capturing code {LARE2D, Arber et al. l200lh . for a simple 2D X-point 
configuration: 



which corresp onds to the magnetic field se en in Figure [T^. Note that Ofman (|l992l '). 
Ofman et al. (| 19931 ) and Steinolfson et al. (| 19951 ) had previously performed nonlinear 
2D calculations of stress ed X-p oints. 

McLaughlin & Hood (|2004l ) and previous results from investigations conducted in a 
cylindrical geometry (!j2]) clearly demonstrate that the Alfven speed (v\ = Bx + By = 
x'^ + y'^ = r^) plays a vital role. Hence, it is natural to consider either a polar coordinate 
system or to drive a circular pulse. In addition, as commented by McClements et al. 
^2004 ). a disturbance initially consisting of a plane wave is refracted as it approaches 
the null in such a way that it becomes more azimuthally-symmetric. Thus, as a first 
step in investigating the nonlinear regime, it is appropriate to consider an azimuthally- 
symmetric initial condition in velocity of the form: 

vj^(2:,y,t = 0) = 2Csin[7r(r-4.5)] for 4.5 < r < 5.5 , (25) 
V|| (x,y,t = 0) ^0 , 

which corresponds to a circular, sinusoidal wave pulse in vj^, where 2(7 is the initial 
amplitude. When the simulation begins, this initial pulse naturally splits into two 
waves, each of amplitude C: an outgoing wave and an incoming wave. The authors 
focus on the incoming wave, i.e. the wave travelling towards the null point, and set 
C = 1. 

The authors find that at early times the incoming wave (identified as a linear fast 
wave) propagates across the magnetic fieldlines and that the initial pulse profile (an 
annulus) contracts as the wave approaches the null point. This is the same refraction 
behaviour that has previously been reported. The authors also note that the incoming 
wave pulse develops an asymmetry, where in the y~ / x— direction the wave peak / 
trailing footpoint is catching up with the leading footpoint / wave peak, and eventually 
forms discontinuities. This can be seen in Figure [5^. The asymmetry develops directly 
due to the choice of a velocity initial condition. In the nonlinear regime, specifying 
an initial condition in velocity also prescribes a background velocity profile. Thus, the 
initial condition (equation I25p appears to excite the m = mode, but this actually 
corresponds to the m — 2 mode in cartesian components. 

At later times, these discontinuities develop into fast oblique magnetic shock waves, 
leading to local heating of the plasma. In addition, the shocks above and below y — 
began to overlap, forming a triangular 'cusp' (called the shock-cusp) and this leads to 
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time 

Fig. 5 Contours of for a fast wave pulse initially located at a radius r = 5 and its resultant 
propagation at (a) time t = I and (b) time t = 2.6. The black lines denote the (changing) 
separatrices and the null point is located at their intersection (origin). Note in (b) that the 
separatrices have been deformed and now form a 'cusp-like' magnetic field structure. The 
amplitude of vx varies substantially throughout the evolution, and hence each subfigure is 
assigned its own colour bar. (c) Plot of time evolution of Jz(0, 0) for < t < 60 . Insert shows 
time evolution of jz(0, 0) for 25 < t < 60 (i.e. same horizontal axis, different vertical axis). 
Dashed lines indicate maxima (red) and minima (blue). Green line shows limiting value of 
jz{0,0) = 0.8615. 



the development of hot jets, which again heat the local plasma and significantl y bent 
the local magnetic fieldlines. The hot jets (which fit the description of Forbes Il988l ') 
set up slow oblique magnetic shock waves emanating from the shock-cusp. In addition, 
there is evidence of slow shocks along the sides of the jet upstream of the tip and 
we see kinks in the fieldlines at the tip of the jet, indicative of a fast shock. Thus, 
the jet heating itself is accomplished by a combination of slow and fast shocks. It is 
interesting to note that the jet has a bimodal structure consisting of a hot, narrow jet 
incased within a broader, lower temperature jet, which is a feature that is not predicted 
by steady-state reconnection theory. 

Eventually, the fast shocks reach the null point and the shocks have deformed the 
magnetic field such that the separatrices now touch one anot her ra ther than intersecting 
at a non-zero angle (called 'cusp- hke' by Priest & Cowlev fl975l ') . This can be seen in 
Figure[S}3. However, the separatrices continue to evolve and so this field structure is not 
sustained for any length of time. The (deformed) null point itself continues to collapse 
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and forms a horizontal current sheet (again, such behaviour was not seen in the linear 
systems). 

Subsequently, the system evolves as follows: the hot jets to the left and right of the 
null continue to heat the plasma, which in turn expands. This expansion squashes and 
shortens the horizontal current sheet, forcing the separatrices apart. The (squashed) 
horizontal current sheet then returns to a 'cusp-like' null point which, due to the 
continuing expansion from the heated plasma, in turn forms a vertical current sheet. 
The evolution then proceeds through a series of horizontal and vertical current sheets 
and displays oscilla tory behaviour, in a similar manner to the linear results of Craig & 
McClymont ( 1993)- The oscillatory nature of the system can be clearly seen from the 
time evolution of jz(0,0) shown in Figure [5j;. The red / blue lines indicate maxima / 
minima in the system and the green line shows jz(0,0) — 0.8615, which is the limiting 
value of the oscillation. 

It is interesting to note that the final state is non-potential, though still in force- 
balance, since there is a finite amount of current left in the system. This is because the 
plasma pressure in the final state is greater to the left and right of the null point than 
that above and below, due to the asymmetric heating from the hot jets. Of course, the 
system will eventually return to a potential state due to diffusion, but this will occur on 
a far greater timescale than that considered by the simulation (idiffusion ~ 

McLaughlin et al l|2009l ') provide two pieces of evidence for reconnection in their 
system. Qualitatively, they observe changes in fieldline connectivity and quantitatively 
they look at the evolution of the vector potential at the null point. Since they have both 
oscillatory behaviour and evidence for reconnection, they conclude that the system 
displays oscillatory reconnection (as detailed by Craig & McClymont Il99ll '). 

The authors then extended the study to look at the effect of changing the ampli- 
tude, C, of their initial condition (equation I25|l . They conclude that a larger initial 
amplitude results in a larger amount of current being left in the system at the end of 
the simulation, i.e. in the non-potential final state. 

Thus, it is clear that the nonlinear behaviour is completely different to that of the 
linear regime. For example, current density now accumulates at many locations, such 
as along horizontal or vertical current sheets, along slow oblique magnetic shocks and 
at the location of shock-cusps. This was not the case in the linear regime, where all 
the current density accumulated at t he nu ll point exponentially in time. 

The work of McLaughlin et al. (|2009h provides a link between two traditionally 
separate areas of solar physics: MHD wave theory and reconnection, and is the first 
demonstration of reconnection naturally driven by MHD wave propagation. 



3.7 Quasi-Periodic Pulsations 

Quasi-Periodic Pulsations (QPPs) ha ve bee n observed in radio, op tical a nd X-ray emis- 
sion of s olar fl ares (e.g. Inglis et al. l2008l : Inglis & Nakari akov 20091: Nakariakov & 
Melnikov IgOOOh and steUar flares (e.g. Mathioudakis et al. l2003l : bood : Mitra-Kraev 



et al. bOOa T QPPs are observed in emission intensity as oscillations with a character- 
istic period (from a few seconds to several minutes) and occur sometimes stably for 
several minutes and at other times in short bursts. Several mechanisms have been pro- 
posed for the ex plana tion of QPPs, either linking the observations to MHD oscillations 
(see Nakariakov 1200?] for a recent review ) or b eing associate d with periodic regimes of 
magnetic reconnection (e.g. Kliem et al. I2OO0I : Ofman & Sui bood ). 
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Fast magnetoacoustic wave behaviour in the neighbourhood ofmagnetic null points 
provides two alternative mechanisms. Firstly, Nakariakov et al. (l2006l l extend the model 
of McLaughlin & Hood 
effect wraps the fast wave around the null point, and large gradients and currents de- 
velop. By utilising a harmonic driver, the resultant current growth and accumulation is 
itself periodic, i.e. the periodicity of the incoming fast waves is efficiently transmitted 
into the periodic modulation of the current density. In addition, Nakariakov et al. pro- 
vide an explanation of where such an incoming fast wave could originate: fast waves can 
leak from an oscillating loop situated near a null point but magnetically disconnected 
from it. Interestingly, this would imply that the period of the observed QPPs is linked 
to the properties of the neighbouring, oscillating loop. 

An alternative explanation may be provided by the work of McLaughlin et al. 
( 20091 ) in which oscillatory reconnection is naturally driven by an incoming fast mag- 
netoacoustic wave. Here, the oscillatory behaviour is due to a cycle of horizontal and 
verti cal current sheets (i.e. a different physical mechanism to that of Nakariakov et al. 
I2OO6I I. However, further work is needed to develop and generalise the model into a di- 
agnostic tool, and such developments should, for example, consider the consequences of 
different initial conditions such, as utilising a localised increase in pressure or internal 
energy in a /3 / plasma, as well as conduct a detailed parametric study. 



20041 ) to include a harmonic driver. As before, the refraction 



3.8 Phase-Mixing 

Finally, we consider a novel piece of work by Fruit fc Craig (|2006l ') whereby the phe- 
nomenon of phase-mixing (Heyvaerts & Priest Il983h is invoked in the viscous and re- 
sistive dissipation of standing Alfven waves within a line-tied X-point geometry. Here, 
the authors consider magnetic fieldlines anchored into a rigid, reflective boundary and 
excite the system with an initial velocity proflle. During the simulation, each fieldline, 
being rigidly tied at the end, oscillates back and forth, with the oscillation frequency 
depending upon the length and magnetic field strength of that fieldline. The phase dif- 
ference between neighbouring lines increases as time evolves, leading to growing cross- 
gradients, and eventually leading to strong visco-resistive damping. Thus, the authors 
conclude that phase-mixing can provide an efficient mechanism for energy dissipation 

of standing Alfven waves in th e vicinity of 2D null points. 

Craig & Litvinenko (|2007l 'l extend the model of Fruit & Craig (|2006l ') to include 



anisotropic (Braginskii bulk) viscosity, and find that the main results are still valid. 



4 Weak guide-field 

An obvious and essential extension of the 2D work described so far is to extend the 
models to full 3D simulations. However, two extension are possible that could provide 
useful signposts before the move to 3D. Firstly, it is possible to extend the model 
to 2.5D with the addition of a third spatial coordinate, by taking into account an 
extra Fourier component of the form e*™^, where m is the azimuthal mode number. 
Secondly, we can consider the addition of a small axial magnetic field perpendicular 
to the plane of the magnetic X-point, i.e. a longitudinal guide-field. Note that with 
the addition of a longitudinal guide-field we are no longer considering null points; even 
though the X-point geometry remains, B 7^ at the X-point (due to the cross field) 
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and we actually have an X-line. However, since the two systems are so closely related 
we will still provide a literature review of the relevant papers. 

Extending the model to 2.5D and 3D will lead to coupling of all the wave modes, 
and thus will most likely result in energy accumulating at both the separatrices and 
the null points. 

The first work in this area was done by Bulanov et al. Bulanov et al. 

considered an X-point magnetic field configuration with a longitudinal (along the X- 
line) magnetic field _B|| . The authors write down, for the first time, the governing 
linearised equations in the ideal /3 = regime, and show that the fast magnetoacoustic 
wave and Alfven wave are linearly coupled by the gradients in the field. Furthermore, 
the authors show that magnetoacoustic perturbations can transform into Alfven waves 
and vice versa, thus leading to current density accumulation at either the null point or 
along the separatrices. However, the rate or efficiency of this process is not explored 
in this analytical work. In the limit of By — >■ 0, the two modes are decoupled and the 
results of 2D work are recovered. 

McClements et al. (|200d ) also investigated the coupling of Alfven waves and fast 
waves in the vicinity of a magnetic X-point with a wea k long itudinal guide field present 
(-By <^ -B_l), and extended the model of Bulanov et al. (1992) to include resistive effects. 
These authors solve the initial value problem for a fast wave being driven by a harmonic 
Alfven wave train and find that energy is channelled into the fast wave, and that large 
gradients start to build-up near the X-point, due to the Alfven speed profile. The 
results indicate that a significant fractio n of the Alfven wave energy is converte d into 
fast wave energy. Ben Ayed et al. (|2009l 'l extend the work of McClements et al. (|2006l ') 
to include a strong guide-field ^ B_l). Again, the authors found that the Alfven 
wave is coupled into the fast mode, with the coupling strongest on the separatrices and 
far from the X-line. 

The three wor ks (Bulanov et al. 1 19921 : McClements et aZ. |200g; Ben Ayed et al. 
l2009|l all assume linear, P — plasma, and McClements et al. (|2006l l and Ben Ayed et 
al. \200di ] concentrate on driving linear Alfven wave trains and observing the coupling 
to the fast mode. 

Conversely, Landi et al. (|2005l ) considered the nonlinear propagation of a harmonic 
Alfven wave train in a 2.5D geometry, consisting of a magnetic X-point threaded by an 
axial magnetic field. Interestingly, the authors find that the driven Alfven waves couple 
to the fast mode through the magnetic geometry, and that the generated fast waves 
have a frequency equal to that of the driven Alfven waves and an amplitude that scales 
linearly with the amplitude of the incoming Alfven waves. This is different to nonlinear 
formation of fast waves from a propagating Alfven wave (due to the ponderomotive 
force), in which the generated fast w aves have a freq uency twice that of the driven 
Alfven wave (Nakariakov et al. Il997l : Tsiklauri et al. I2OO11 : 120021 ') and an amplitude 
related to the square of the driven Alfven wave amplitude. Thus, Landi et al. ( 20051 ') 
propose that this indicates a mechanism of mode conversion that differs from the 
standard nonlinear fast wave excitation via the ponderomotive force. The authors also 
find that these generated fast waves rapidly develop into fast-mode shocks, and thus 
the wave dissipation is concentrated into thin current structures. However, this is not 
surprising as the authors consider very large amplitudes for their driven Alfven waves, 
precisely so that they can observe shock formation within their simulation domain. 
The authors do not report whether the fast waves experience the refraction effect or 
whether the Alfven waves are confined to the magnetic fieldlines that originate on. 
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McClements et al. (|2006t ) and Be n Aye d et al. (|2009l l consider radial symmetry for 
their driven waves, and Landi et al. (|2005h considers a large-amplitude, periodic wave 
train driven simultaneously on both side boundaries. Hence, it would be interesting to 
see if the results of all these papers persist if one considers a single wave pulse, as this 
would allow the transient behaviour to be observed, and also if the waves were driven 
on a single side of the numerical domain. 



5 Three-dimensional null points 

Finally, let us now review the behaviour of MHD waves in the neighbourhood of 3D 
null points. However, surprisingly few papers have been written that address this issue, 
or at least, papers that concentrate on the the transient propagation of the modes. 
Most papers have focused on the dynamics of current formation in an attempt to 
locate regions where reconnection is most likely to occur, rather than on the transient 
propagation of the MHD waves. 

The first study of the d ynam ics of current formation at 3D null points was per- 
formed by Rickard & Titov (|l996l ) . These authors solved the linear, /3 = MHD equa- 
tions, and studied a 3D null point that is axisymmetric about the spine (i.e. a proper 
3D null, see i31.3.2|l . Perturbations were decomposed into azimuthal modes labelled by 
mode number m, and the analysis is performed in cylindrical geometry. 

Several different perturbations are considered by Rickard & Titov (see their figure 
2) and are driven with a velocity pulse (in Vr,vg,Vz) on the boundaries (either the radial 
boundaries or upper/lower boundaries). Numerical simulations show that axisymmetric 
perturbations (m = 0) lead to current accumulation along the spine, whilst the m = 1 
mode produces currents in the fan plane and at the null point itself. For ?n > 1, 
there is no current accumulation anywhere along the skeleton. Rickard & Titov also 
noted that in these axisymmetric equilibria, the azimuthal components decouple from 
the remaining components when m — 0, and that in 2D only the m = mode was 
associated with producing currents at the null. 

The primary aim of this study was to investigate current accumulation, and not to 
investigate the behaviour of MHD waves around 3D nulls point. Of course, the pertur- 
bations investigated correspond to a combination of Alfven wave and fast magnetoa- 
coustic waves, but the authors do not use this terminology. In addition, the authors 
drive a combination of Vr,vg,Vz which generates both wave types, and thus in some 
cases it is difficult to confirm that current accumulation results from a certain wave- 
type. However, some behaviours are clear: a pure m = wave pulse corresponds to 
a torsional Alfven wave, which the authors note is channelled along the equilibrium 
magnetic fieldlines. Secondly, the authors note that the m = 1 motions can propagate 
across magnetic fieldlines, and that a focusing effect is seen in the evolution of jr and 
jg (as we would expect for the fa st wave) for certain disturbances. 

Thus, Rickard & Titov l|l996h give a tantalising suggestion that our understanding 
of 2D behaviour of the fast and Alfven wave transfers to 3D. However, it is not possible 
to see from their paper if driving a pure Alfven wave results in generation of a fast 
mode disturbance through the magnetic geometry, or vice versa, as expected from ^31 
The simulations of Rickard & Titov also utilise refiecting boundary conditions in their 
cylindrical geometry. In addition, it can be argued that their investiga tion i s actually 
2.5D and not a fully 3D experiment. Nevertheless, Rickard & Titov ( 19961 ) is still a 



landmark paper for the investigation of 3D MHD wave behaviour in the neighbourhood 
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ofjnagnetic null points. It would be interesting to repeat the work of Rickard & Titov 
(|l996l ) but to try to excite pure modes and to focus on the resultant transient wave 
propagation and possible mode conversion. 

The work of Rickard & Titov (| 19961 ') has been extended to multiple null point 
topolo gies in a series of papers by Galsgaard and co-workers. Galsgaard et al. (| 19961 : 
Il997al l looked at shearing a 3D potential null point pair, with continuous (opposite) 
shear on two opposite boundaries (parallel to the separator), where the fieldlines were 
not returned to their original position. This generated a wave pulse that travelled 
towards the interior of the domain from both directions, and resulted in current accu- 
mul ation a long the separator line with maximum value at the null points. Galsgaard et 
al. (|l997bh looked at perturbations in 3D magnetic configurations containing a double 
null point pair connected by a separator. The boundary motions used were very similar 
to those described above (i.e. shear the boundary and fix). Their experiments showed 
that the nulls can either accumulate current indivi dually or act together to produce 
current along the separator. Galsgaard & Nordlund ( 19971 ) found that when a magnetic 
structure containing eight null points is perturbed, current density accumulates along 
separator lines. 

In all these Galsgaard et al. papers, the boundary conditions tried to mimic the 
effect of photospheric fo otpoin t motions by moving the boundary and holding it fixed. 
As in Rickard & Titov |1996|), the terminology of MHD waves was not invoked, and 
the perturbations considered were a combination of MHD waves. 

The first investigation specifically looking at MHD wave propagation about a proper 
3D null point was reported in Galsgaard et al. \200^ ). where the authors looked at a 
particular type of wave disturbance and solved the nonlinear /3 = MHD equations. 
Galsgaard et al. investigated the effect of rotating the field lines around the 

spine to generate a twist wave (essentially a torsional Alfven wave) and followed its 
propagation towards the null point. Twists were imposed simultaneously on the upper 
and lower boundaries, with both the same and opposite vorticities considered. The 
authors found that the helical Alfven wave spreads out as it propagates towards the null 
point and is confined to the magnetic fieldlines it originates on. As the wave approaches 
the fan plane, the wave spreads out along the diverging fieldlines and produces current 
accumulation in the fan plane. 

In addition, the authors also observe the generation of a fast-mode wave, and find 
that this wave focuses and wraps around the null point. The authors suggest that 
their twist wave is a pure helical Alfven wave, and that nonlinear effects generate the 
fast-m ode wave (using the ponderomotive mechanism detailed by Nakariakov et al. 
Il997l l. However, the authors themselves also note that the fast wave appears to be 
absent where the boundary driving is slowly increased and only appears when a near- 
discontinuity in the boundary driving velocity is utilised. Thus, the exact nature of the 
fast-mode generation is unclear. 

Galsgaard et al. confirm that several of the key properties of fast and Alfven 

waves transfer to 3D geometry. However, again their investigation is actually 2.5D and 
not fully 3D. In addition, closed boundaries were used, and thus waves propagate 
toward the boundaries and are refiected back into the domain. 

The authors also analyse the linear /? = MHD equations using the WKB method 
for their azimuthally-symmetric 3D null, and find good agreement. They find that 
the equations for fast and Alfven perturbations decouple, although this is not sur- 
prising as, due to their choice of symmetrical geometry, their resultant equations are 
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two-dimensional (since a proper 3D null point is essentially 2D in cylindrical polar 
coordinates). 

Pontin & Galsgaard (|2007h and Pontin et al. have performed numerical 

simulations in which the spine and fan of a proper 3D null point are subjected to 
rotational and shear perturbations. They found that rotations of the fan plane lead to 
current density accumulation about the spine, and rotations about the spine lead to 
current sheets in the fan plane. In addition, shearing perturbations lead to 3D localised 
current sheets focused at the null point itself. Again, this is in good agreement with 
what we may expect for MHD wave behaviour, i.e. current accumulation at certain 
parts of the topology. 

The first study of MHD waves in t he nei ghbourhood of an improper 3D null point 
was investigated by McLaughlin et al. (|2008l ). These authors consider the propagation 
of the fast and Alfven waves about the magnetic equilibrium Bq = [x, ey, — (e -f 1) z], 
for both proper (e = 1) and improper (e > 0, e 7^ 1) 3D potential null points. The 
authors utilise the 3D WKB approximation of the linear /? = MHD equations. 



McLaughlin et al. (|2008l ^ find that the fast magnetoacoustic wave experiences re- 
fraction towards the magnetic null point, and confirms that, as in 2D, the effect of this 
refraction is dictated by the Alfven speed profile. The fast wave, and thus the wave 
energy, accumulates at the null point. The current build-up is shown to be exponential 
and the value of the exponent depends upon e. Thus, as in 2D, there is preferential 
heating at the null point for the fast wave. 

For the Alfven wave, the authors find that the wave propagates along the equilib- 
rium fieldlines and that a fluid element is confined to the fieldline it starts on. For an 
Alfven wave generated along the fan-plane, the wave accumulates along the spine. For 
an Alfven wave generated across the spine, the value of e determines where the wave 
accumulation will occur: either the fan-plane (e = 1), along the a;— axis (0 < e < 1) 
or along the j/— axis (e > 1). Analytical results show that current density builds up 
exponentially, leading to preferential heating in these areas. 

McLaughlin et al. (|2008h also provide a quantitative analysis of the preferential 
heating/energy release process. They derive an analytical expression for current evo- 
lution resulting from fast wave propagation along the spine (their equation 23) which 
is of the form |j| = oje^'^^^-'*/ [zq (e -|- 1)] (where ui is the wave frequency and zq is the 
starting point on the spine). Furthermore, in order to give an order-of-magnitude esti- 
mate, the authors show that for characteristic coronal values (L = 10 Mm, i? = 10 G, 
po ~ 10^^^ kg m"'^) a planar fast wave propagating along the spine will build-up a 
current of 0.3 mA after a time of t = 1 seconds, and that resistive effects become non- 
negligible after a time t = log (cjA^/ry)/4 ~ 7 seconds (where after 7 seconds, the fast 
wave has built up a current of 0.87 mA and has travelled a distance of 7.13 Mm). The 
Alfven wave is degenerate with the fast wave along the spine in a /? = plasma, and 
so has identical estimates for these characteristic conditions. 

However, McLaughlin et al. (|2008l ) are unable to reach any conclusions about the 
coupling of the fast and Alfven wave types due to the geometry of the magnetic field. 
The WKB approximation, in the form utilised by McLaughlin et al., does not take into 
account the coupling of the fast and Alfven wave types due to the geometry of the 
magnetic field. Under their approximation, the waves actually see the magnetic field 
as locally uniform. 

Thus, a fully 3D model of MHD wave behaviour in the neighbourhood of a general 
3D null point has yet to be investigated, where the investigation tracks the propagation 
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and evolution of each MHD wave, and determines the efficiency of mode-conversion due 
to the magnetic geometry against that due to nonlinear effects. 

6 Conclusions 

The behaviour of all three MHD wave types; Alfven, fast and slow wave, has been 
investigated in the neighbourhood of 2D, 2.5D and (to a certain extent) 3D magnetic 
null points, in a variety of geometries and under a variety of assumptions. The main 
conclusions may be summarised as follows: 

— The linear, fast magnetoacoustic wave behaviour is dictated by the equilibrium fast 
wave speed profile (i.e. y^u^ + Cg), which in low-/3 plasmas can be thought of as 
the equilibrium Alfven-speed profile. The fast wave is guided towards the null point 
by a refraction effect and wraps around it. The fast wave slows as it approaches the 
null, leading to a decrease in length scales and thus an increase in current density 
close to the null point. In a /? = plasma, the fast wave cannot cross the null point 
and the build-up of current is exponential, indicating that dissipation will occur on 
a timescale related to log?7. Thus, linear fast wave dissipation is very efficient, and 
null points mil be locations for preferential heating. For /3 7^ 0, the fast wave can 
cross the null point, due to the finite sound speed there, and wave energy can now 
escape the null point. In this case, there exists two competing phenomena and the 
dominate effect is determined by the value of the plasma- /3. 

— The linear Alfven wave propagates along the equilibrium fieldlines and a fiuid el- 
ement is confined to the fieldline it starts on. Since the propagation follows the 
fieldlines, the Alfven wave spreads out as it approaches the diverging null point. 
In 2D, all the Alfven wave energy accumulates along the separatrices and the cur- 
rent build-up is exponential in time. In 3D, for an Alfven wave generated along 
the fan-plane, the wave accumulates along the spine and for an Alfven wave gen- 
erated across the spine, the value of e determines where the wave accumulation 
will occur: fan-plane {e — 1), along the 2;— axis (0 < e < 1) or along the y— axis 
(e > 1). Hence, all the Alfven wave energy will be dissipated along the separatrices/ 
separatnx surfaces and these will be the locations for preferential heating. 

— The behaviour of the slow wave in the neighbourhood of null points has received the 
least attention in the literature. The linear slow wave is found to be wave-guided 
and accumulates along the separatrices. A low-/3 fast wave can generate/convert 
into both a high-/? fast and high-/? slow wave as it crosses the = eg mode- 
conversion layer. Such a layer is a natural consequence for a null point emersed in 
a /3 plasma. In fact, the value of /3 grows as close to a null point. 

— The addition of a weak guiding field leads to linear coupling between the fast and 
Alfven waves in a low-/3 plasma, and thus the propagation of either mode can 
generate the other. Such a configuration is, of course, no longer a null point, but 
rather an X-line. However, the nature of mode-coupling for 3D null points is, at this 
time, uncertain. Fast waves have been shown to be generated by the propagation of 
the Alfven wave, but it is unclear if this is due to the fieldline geometry, nonlineax 
coupling or both waves being simultaneous generated by a common driver. 

— Results in 2D show that in the nonlinear regime, the fast magnetoacoustic wave can 
deform the equilibrium X-point configuration, leading to a cycle of horizontal and 
vertical current sheets and associated changes in connectivity. Thus, the system 
exhibits oscillatory reconnection. 
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It is clear that the equilibrium magnetic field plays a fundamental role in the propa- 
gation and properties of MHD waves. In general, an arbitrary disturbance/perturbation 
will generate all three wave modes and current accumulation could occur at all the 
null points, and/or along the spine, fan and separators. Thus, the results described 
in this review all highlight the importance of understanding the magnetic topology 
in determining the locations of wave heating. 
However, several big questions still remain in this area: 

The nature of the coupling of the three modes in 3D needs to be addressed, and 
the importance of coupling due to the magnetic geometry verses nonlinear coupling 
should be investigated. 

The theory of nonlinear fast waves driving oscillatory reconnection should be ex- 
tended to study more general distur bances, and to investigate how robust the initial 
findings of McLaughUn et al. (|2009l ') 



are. 



— The key results for the linear fast and Alfven wave make clear predictions as to 
where preferential heating can occur. It would be interesting to see the t heore tical 
models developed wit h forw ard modelling (see e.g. Kilmchuk & Cargill I2OO1I : De 
Moortel & Bradshaw bOOSi ) to provide tell-tale observational signatures, and for 
these synthetic results to be compared with observational data. 

In conclusion, we have seen that the study of MHD wave behaviour in the neigh- 
bourhood of magnetic null points is a fundamental plasma process, and can provide crit- 
ical insights into other areas of plasma behaviour including: mode-conversion (i ]3.5|) . os- 
cillatory reconnection fi )3.6|) . quasi-periodic pulsations ( i)3.7|) and phase- mixing f i)3.8|) . 

W e now know that the corona is full of MHD wave perturbations (Tomczyk et al. 
We also know that null points are an inevitable consequence of the distributed 
isolated magnetic fiux sources at the photospheric surface, and potential and non- 
potential field extrapolations suggest that there are always likely to be null points in 
the corona (see t|1.4|l . Thus, these two areas of scientific study (MHD wave behaviour 
and magnetic topology) will inevitably encounter each other at some point, i.e. MHD 
waves will propagate in the neighbourhood of coronal null points. Thus, MHD wave 
propagation about magnetic null points is itself - theoretically - a fundamental coronal 
process. 

However, there is as yet no clear observational evidence for MHD wave behaviour 
in the neighbourhood of coronal null points. In the lead author's opinion, the successful 
detection of MHD oscillations around coronal null points will require input from two 
areas: high-spatial/temporal resolution imaging data as well as potential/non-potential 
extrapolations from co-temporal magnetograms. Two of the instruments onboard the 
recently launched Solar Dynamics Observatory (SDO) may satisfy these requirements: 
the Atmospheric Imaging Assem,bly (which will provide high-quality imaging data) 
and the Helioseismic and Magnetic Imager (which will provide vector magnetograms). 
Thus, the first detection of MHD waves in the neighbourhood of coronal null points 
may be reported in the near future. 
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